Spin-orbit interaction from low-symmetry localized defects in 
semiconductors 
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Abstract -The presence of low-symmetry impurities or defect complexes in the zinc-blende direct- 
gap semiconductors (e.g. interstitials, DX-centers) results in a novel spin-orbit term in the effective 
Hamiltonian for the conduction band. The new extrinsic spin-orbit interaction is proportional to 
the matrix element of the defect potential between the conduction and the valence bands. Because 
this interaction arises already in the first order of the expansion of the effective Ifamiltonian in 
powers of Uaxt/Eg <C 1 (where f/cxt is the pseudopotential of an interstitial atom, and Eg is the 
band gap), its contribution to the spin relaxation rate may exceed that of the previously studied 
extrinsic contributions, even for moderate concentrations of impurities. 



Introduction. — The spin-orbit interaction (SOI) is 
a remarkable example of an "emergent interaction" . For 
electrons in vacuum, the SOI emerges from the projection 
of the Dirac equation on positive energy states. It is in- 
trinsically weak, since its "coupling constant" contains the 
electron rest energy mc^ in the denominator, and thus re- 
quires strong nuclear potentials to become truly effective. 
For electrons in the conduction band of semiconductors, 
the situation is different. The spin and orbital parts of the 
Bloch wave functions for these electrons are already entan- 
gled by the spin-orbit interaction arising from the strong 
periodic potential in the atomic cores. Therefore, the po- 
tential from a symmetry-breaking defect leads to a SOI 
that is orders of magnitude larger [TH3] than the SOI one 
would expect in vacuum for the same strength of the po- 
tential. This enhanced SOI plays a crucial role in recently 
proposed schemes for the electric manipulation of spins 
in semiconductors [4l[5]. So far, however, the attention 
has been focused on defects that break the translational 
symmetry of the lattice (e.g. substitutional defects, va- 
cancies) leaving the point group intact. But the breaking 
of translational symmetry, as will be shown below, adds 
the defect potential C/oxt to the diagonal elements of the 
Hamiltonian [21 IS], whose conduction and valence states 
are already split by the band gap Eg 3> t/cxt, and leads 
only to relatively small corrections (the so-called "central 
cell corrections" [7]) to energy spectra, effective electron 



mass, etc. Contrary to that, the point-group symmetry re- 
duction results in a correction to the off- diagonal matrix 
elements (between the conduction and the valence band) 
of the Hamiltonian, which are small to begin with; conse- 
quently the corresponding corrections are not expected to 
be relatively small, especially in a wide-gap semiconduc- 
tor where one might expect the conduction-band SOI to 
be more sensitive to a lowering of the point group symme- 
try of an atomic cell than to the breaking of translational 
symmetry. 

In this Letter we demonstrate that lowering the point 
group symmetry leads to a novel type of spin- orbit interac- 
tion [see Eq. (|14p below] in direct-gap semiconductors, and 
we explore some implications for the practically important 
calculation of the spin relaxation rate plIMIS] . The reduc- 
tion in point group symmetry provides the leading-order 
contribution to spin relaxation whenever the conventional 
D'yakonov-Perel' (DP) mechanism [6] due to the spin-split 
band structure is absent (this happens, for example, in 
semiconductor quantum wells grown along certain direc- 
tions [H]). We will demonstrate this point by explicitly 
comparing the spin relaxation rate produced by interstitial 
defects with the spin relaxation rate produced by substitu- 
tional defects with a comparable elastic relaxation time r. 
Furthermore, even though the DP mechanism is frequently 
thought to be the dominant mechanism for spin relaxation 
in III-V semiconductors (both in the bulk and in quantum 
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wells) pi?lll2lll3| , we find that in some cases the dominance 
of DP is not so clear-cut: our mechanism can compete and 
even dominate for strong enough disorder, so long as the 
disorder docs not invalidate the diffusive transport regime. 

A direct-gap semiconductor subjected to an extrinsic 
field Ucxt{r) is described by the Hamiltonian 



2m 



(1) 



where p = —ihVr is the momentum operator, Uat is the 
periodic crystal potential and t/ext is an "extrinsic" poten- 
tial from impurities, defects, and/or other sources. 

In the standard k -p approach, which we adopt here, the 
Hamiltonian is represented as a matrix on the basis of 
wave functions i'^^ir) = e^'^'^Unoi'r), where Uno{f) are 
the periodic parts of the Bloch wave functions of the peri- 
odic crystal at fc = 0. Notice that ip^j: are not Bloch wave 

functions when k ^ 0: for this reason, even in the absence 
of the external potential, the Hamiltonian has off-diagonal 
matrix elements of the form k ■ Pnn'5j:p /m between ■0^^ 
and V'n'feo with n ^ n' . Here p„„/ is the matrix element 
of p between Un^(f) and u,i'o(?^). The field UaAif) lowers 
the otherwise perfect symmetry of the crystal in two ways: 

1. It breaks the translational symmetry by creating 
diagonal (intra-band) matrix elements of the form 
Ucxt{k — k')Snn', where Uextik) is the Fourier trans- 
form of Ucxtir); 

2. It may also lower the local point group symmetry, 
producing additional off-diagonal (inter-band) matrix 
elements with n ^ n' and k = k' . 

An interstitial atom or a point defect such as a DX-center 
[7] or an EL2-complex [15] affects both symmetries: the 
translational one and the point-group one. The traditional 
approach to the generation of spin-orbit interaction fo- 
cuses on the intra-band part of Ue^t, which arises from 
the breaking of translational symmetry. Here we concen- 
trate on the inter-band part of C/cxt , which arises from the 
lowering of point-group symmetry. 

Let us consider a situation in which the extrinsic 
field Uext is produced by short-range interstitial defects, 
e.g., DX-centers, which are known to appear sponta- 
neously in AlxGai_xAs for 0.22 < x < 0.5 [TB]. We assume 
that the field induced by every single defect is negligibly 
small outside of the atomic cell in which the defect resides. 

This assumption makes it convenient to separate the 
coordinate variable r (determining an arbitrary point 
inside the crystal) into a "macroscopic" (p) and a 
"microscopic" (r) part: p is a discrete variable pointing 
to the origin of each unit cell, while r varies within the 
cell. Any point inside the crystal is thus uniquely iden- 
tified by the pair of variables {p,f) so that f ~ p + f. 
With this notation our external potential can now be fac- 
torized into the product of a macroscopic part, U{p), and 



a (dimensionless) microscopic part, V{r): 

u,xt{p,f) = u{p)v{f), f y(f)dV = \vi (2) 

where v denotes the unit cell centered at the origin, and \v\ 
is its volume. Let us also assume that the defects are 
placed into a randomly chosen set of unit cells, centered 
at positions pj with j ^ 1 . . . n. Then we can write 



C/(P) = C/(0)^X(P-P,), Xip] 



i/\v\, pev, 

0, p^v, 



(3) 



where n is the total number of defects, and C/*^*^^ is the 
amplitude of the potential. 

In the macroscopic limit (which corresponds to the elec- 
tron gas model) \v\ — >■ 0, so we have xiP~ Pj) ~^ ^iP~ Pj) 
in ^ and the disorder potential is described by the stan- 
dard Lifshitz disorder model |17| 



Uip)U{p') = 



5{P- P') , 



(4) 



where ve is the density of states at the considered electron 
energy, te is the elastic momentum relaxation time at the 
same energy, and the overline stands for the averaging over 
different disorder configurations. 

In what follows we will be concerned with (direct-gap) 
zinc-blende semiconductors, for which the two-fold degen- 
erate conduction band is described by the Fg-irreducible 
representation (irrep) of the double cubic group, while the 
six-fold degenerate valence band is a mixture of two sub- 
bands described by irreps Fg and Fy. 

In this Letter we use the 8 x 8 • p model, where the 
initial Hamiltonian ([!]) is reduced to the 8 x 8-matrix Hp 
whose elements are operators acting on p-dependent func- 
tions (where pis a long-range variable defined above). Let 
us divide the 8x8 matrix Hp into four blocks as follows: 

m 



Hy 




(5) 



He Hcv 

He ~ HyS + Eg, HyS ~ H yj + A, 

where A is the spin-orbit-induced splitting between the 
Fg and F7 valence bands; the blocks He and Hyr are two- 
dimensional, and the block Hyg is four-dimensional. Using 
the well-known mathematical formula for the inverse of a 
block matrix, we express the 8x8 Green's function (GF) 
= [E — H]^^ in terms of the matrix blocks defined 
in dSl). 

When holes in the valence band arc absent, all transport 
properties are controlled by electrons in the conduction 
band. Under these conditions the relevant information 
from the full GF is contained entirely in the upper- 
left 2x2 CB-block G!? of G^: 



E-He 



He 



E - Hy 



-1 

H^, 



(6) 
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The conduction-band GF Gf corresponds to the effective 
conduction-band Hamiltonian according to the usual rule 



Gf = [E-H 



=E-{G 



E\-l 



(7) 



The spin-dependent part of the Hamiltonian is en- 
tirely contained within the last term in (|6]). In the ab- 
sence of spin-orbit splitting of the valence bands (i.e., for 

A = 0) would be spin-independent since the following 
sum (over all considered [19j states) is spin-independent: 



i,TlGr7,8 



E — Hr 



\n){n\H, 



(8) 



where the sum is taken over envelope states in the valence 
bands. [We note that due to the diagonal form of in ([S]) 
only terms with / = n survive in (|5]).] Due to the sum 
rule dH) it is possible to express the spin-dependent part 

of as a sum over either r7 or Fg states. We prefer to 
sum over F7 states since this results in simpler expressions; 
then the spin-dependent part of is entirely contained 



E 

ner7 



kp 



+ C/oxt(r') 



where (G^ is a scalar GF: 



kp 



+ f/oxt(r) 



(9) 



(10) 



and the energy zero corresponds to the bottom of the 
conduction band. Different spin-orbit contributions can 
be obtained by expanding (G^-^s-^ — in © in 

powers of {Tj \U,^,\ Tj) /Eg < 1 and (Fy \U,^t\ T7) /{Eg + 
A) ^ 1, where p-dependent (F7 |f/ext| F7) is a 2 x 2 matrix 
obtained by projecting the field C/ext onto the Fy-states. 
Thus in our calculation the quantity (F7 | Uext \ F7) must 
be understood as the pseudopotential associated with the 
point-group-symmetry-breaking defect; the importance of 
these symmetry-breaking potentials extends to regimes of 
larger C/ext (p), for which the higher-order terms in the ex- 
pansion of Gf can sometimes be analytically summed jl5| . 

The novel spin-orbit interaction. — The zeroth 
order of the expansion corresponds to the approximations 
\E—Ef\ ^ Eg (where Ep is the Fermi energy) and C/cxt ^ 
Eg in (G, or, in other words, 

' ' (11) 



<G 



E-E„-A 



- <G 



E-E„ 



E, 



Ea + A' 



^9 

Substituting pT|) in © and disregarding all spin- 
independent terms one obtains 



^^SOI — 77 



1 



1 



S„ + A 



(12) 



where Hc,v7 is the rightmost 2x2 part of Hcv defined in (O 
(in other words, Hc^v7 is the Hamiltonian matrix taken 
between the conduction band states Fg and the valence 
band states F7). Making use of Eq. ([2]) for the external 
potential these matrix elements work out to be 



[Hc,v7]kP 



(ra\e"''-^He''''-^\r7 



(13) 



where P = P7 ^ ih{S\p^\X7) /m and V = (^S\V\R 
Here 15*) denotes the orbital part of the Fg state, and 
the components of the (spin-independent) vector function 
1^) = (IX7), IF7), IZj)) participate in the basis formation 
for the irrep F7. Lastly, the vector a = (cti, (72, fxa) denotes 
the set of three Pauli matrices. From Eqs. (fT2|) and ([T3| 

we extract the spin-dependent part of H: 



i^SOI = -i-r s 
h 



Vx Vp-C/(p) + C/(p)Vp- 



where s 



h/7/2 is the spin operator and 
2P / 1 1 



h = 



En 



En 



(14) 



(15) 



^9 ^9 

Notice that in order to arrive at Eq. (fT4)) we replaced the 
wave vector k of Eq. by the operator Vp. This is legit- 
imate because the basis functions of our representation de- 
pend on k only through the plane wave factor e'*^ '' which, 
in the macroscopic limit, becomes e'*^ '', since k ■ f <^ 1. 
Thus, on this basis, the operator \I p is equivalent to mul- 
tiplication by ifc. 

The Rashba spin-orbit interaction. The first 
order of the expansion of (G^^^s^'^ — ^^^^n in ([9]) in 
powers of C/oxt corresponds to the approximation 

1 1 



(G 



(G 



E-E, 



1 



^(rrlC/cxtlF, 

(F7|C/cxt|F7) 



A' 



(16) 



where (F7 | [/oxt | F7) is a 2 x 2 matrix obtained by projecting 
the field C/ext onto the F7-states. Similarly to how it has 
been done in the previous section, but using (|16p instead 
of (fTTj). one reproduces the Rashba spin-orbit term: 



2h? 



{ka)U{p){ka) 



2P2 
~3~ 



E, 



{VpU{p)) X k 
-'-{Eg + A)-' 



(17) 



The Rashba result (fT?)) is reproduced in case when C/oxt (r) 
is a smooth function so that in ([2]) V{f) = 1 and C/ext(?^) = 
C/cxt (p, f) = U (p ). In a more complicated situation, when 
C/oxt may significantly change on the scale of a single 
atomic cell, central cell corrections arise, which are, how- 
ever, small compared to both (|14p and (IT71) . 
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Discussion. — The values of the SOI lengths h and I2 
for several direct-gap A3B5 semiconductors are listed in 
Table [TJ The spin-orbit interaction in (|14p is different 
from the previously studied extrinsic SOI (|17p in that (i) 
its amplitude contains only one power of Eg in the de- 
nominator and (ii) it occurs due to defects that violate 
the point-group symmetry of an atomic cell. Point (ii) is 
a consequence of the fact that the matrix elements V, de- 
fined in Eq. ([T3|) are essentially null if the "form factor" 
V{f) of the impurity potential has the symmetry of the 
point group of the lattice, as is the case for a substitu- 
tional impurity in a zincblendc semiconductor. Thus, the 
spin-orbit interaction associated with substitutional im- 
purities is of the conventional Rashba-type p7)) . while the 
spin-orbit interaction associated with symmetry-lowering 
defects has an additional inter-band term. 

In the rest of the article we see how the new spin-orbit 
term (|14p contributes to the spin relaxation. 



Electron spin relaxation. — We assume that the 
electrons are initially spin-polarized in the z-dircction, 
with spin distribution function Je = ^ctq + ^^rs- Fol- 
lowing the approach outlined in Ref. [20] we obtain the 
dynamical equation for f'^^: /|?^ = — /^^^/rsoi(£^), where 
the relaxation rate TgQj^[E) is given by 



X Tr 



dEdio 



dt'e-"^^'/^ 



(18) 



where Gf/^ = 



E + h^Wy2m - U{p) ± iO 



are the 



Green's functions operators for the electron gas in the 
presence of disorder ([2]). We calculate the relaxation 
rate ([T8|) using the disorder averaging technique [17] to- 
gether with the loop expansion |21| . i.e., the expansion 
in powers of H/Ete ^1. In a degenerate electron gas 
fEt^O only in the close vicinity of the Fermi energy Ep, 
so that only electrons with energies E close to Ep partic- 
ipate in the relaxation. In this case, the spin relaxation is 
exponential with the relaxation rate t^Jj(£'f): 



/>oo 

-ft / /E^JyETs^i{E)dE 

Jo 



=:'^SOi(^f)- 



(19) 



Fig. [T] shows the relevant diagrams in the zeroth order 
of the loop expansion |21| . Both diagrams have a com- 
mon part: a Green's function with two spin-orbit vertices 
connected by a disorder-averaging (dashed) line. First we 
calculate the integrals over u and t from this common part: 



i Tr 5m / dt 



duj 
2^' 



-iujt /h ^ 



r-2 



2te 



(20) 



where we denoted xsoi = Vzpli/H <C 1, and p = v2rnE. 
After the integration in (|20l) , the two triangle diagrams in 





(a) 



Fig. 1: Two diagrams for the kinetic term in Eq. psp in the 
leading order of h/^Er) <C 1 expansion (zero loop) [21]. The 
solid lines denote Green's functions, the dashed lines represent 
averaging over the different disorder realizations; the ladder 
composed by the dashed lines stands for the cooperon |17j . 
The solid circles represent the new SOI interaction of Eq. (|14p . 



Fig.[l]turn into two bubbles; from (fTS)) and (|20p wc obtain 



'SOI 



L ^A , 



(21) 



so that the spin relaxation is: 

in 2D rsoVi?) = -^K', 
4mEll 



in 3D 



(22) 



where r is the elastic relaxation time. Here by "2D" we 
mean the two-dimensional electron gas in a quantum well 
with only one subband occupied by doped electrons, while 
"3D" refers to the usual bulk case. 

For quantum wells grown along certain particular direc- 
tions (eg., the [110] direction in GaAs) intrinsic spin pre- 
cession is absent for spins oriented along a specific direc- 
tion, and the extrinsic relaxation mechanism dominates. 
The ensuing spin relaxation rate 2V^ppllT~^ /h"^ is pro- 
portional to the elastic scattering rate = t^^. It is in- 
structive to compare the spin relaxation time arising from 
the novel SOI (fT4)) with what would be obtained from the 
conventional approach [SKH] , in which the spin-orbit inter- 
action is given by 



Hi 



SOI, 2 



Hi 



VpU{p) X Vp- /h, 



(23) 



with I2 



mechanisms give a spin relaxation rate proportional to 
T~^, so the ratio of the two contributions is independent 
of T and given by 



tsoi/tsoi.2 - lt/{2a^ll) ~ 8 x 10" 



(24) 



for a GaAs quantum well of width a = lOOA. We see that 
in 2D the relaxation rate due to the novel SOI is two orders 
of magnitude larger than the "conventional" one. In the 
3D case the ratio tsoi/t"soi,2 becomes energy-dependent: 



TSOI 
TSOI.2 



mElj 



0.04 at i;=0.03cV. 



(25) 
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Finally, let us compare our spin relaxation time with the 
well-known DP relaxation time, which is caused by the 
intrinsic spin splitting of the conduction bands due to bulk 
inversion asymmetry (BIA). Our calculation in this case 
gives 



327^4 

in 2D Tb/a (E) = TE^mEb^ , 
64 b'^m^E^ 



(26) 



in 3D 



^bia(^) 



105 



where the parameter b = bfl^'^ is defined in Eq. (6.2a) 
of Ref. [2] and its value is listed in Table I. Comparing 
Eqs. (22) and we see that the ratio of the BIA and 
SOI relaxation rates is given by 



^bia(^^) 



'SOI 



[CETE/hf 



where the constant C = (V 1536/105) (6m//i/i^) is listed 
in the last line of Table I for a few semiconductors. A 
short elastic scattering time favors our extrinsic mecha- 
nism; on the other hand te cannot be made too small in 
our calculation since our diagrammatic approach is valid 
only for weak disorder, Ete > h. The value of the con- 
stant C is reported in Table I for several bulk zinc-blende 
semiconductors. Since C < 1 in all cases, we can conclude 
that the novel SOI will dominate the spin relaxation 
rate in a window of parameters such that the condition 
1 < ETE/h < C^^ is satisfied. As for the 2DEG in quan- 
tum wells, we find that the DP-BIA mechanism generally 
dominates the spin relaxation, except in those special cases 
in which it is absent for symmetry reasons. 

We now estimate some spin relaxation times for GaAs 
in specific situations. For a 3D donor concentration of 
= 1.5 X 10^^ cm-3 we find (see Fig. 3-23 in Ref. [22]) 
that the corresponding mobility is ^ — 4 x 10^cm^/(Vs), 
which corresponds to a 3D clastic relaxation time r = 
/im*/e = 1.5 X 10^^'^s. If all these donors were interstitials, 
with li given above in Table [1] then for a quantum well 
with a Fermi energy of 34meV, the spin relaxation time 
would be 9 X 10~"s. In bulk GaAs at Ef = 26meV we 
estimate the spin relaxation time from this mechanism 
under these same assumptions as 2 x 10" 
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s. 

Conclusion. — The reduction of the local point- 
group symmetry of an atomic cell caused, e.g., by inter- 
stitial atoms, leads to a new type of spin-orbit interac- 
tion [Eq. ([T4| ] which scales inversely to gap size in direct- 
gap zinc-blende (A3B5) semiconductors. We speculate 
that low-symmetry defects may also play a role in enhanc- 
ing the spin-orbit interaction in group IV semiconductors, 
and at polar interfaces between materials of cubic symme- 
try, such as at the LaAlOs/SrTiOs interface [23]. Other 
quantities, which originate from off-diagonal matrix ele- 
ments of the k • p Hamiltonian {e.g., the renormalization 
of the electron mass and of the g factor) will also be af- 
fected. 





GaAs 


GaSb 


InP 


InAs 


ZnSe 


CdTe 




1.52 


0.73 


1.34 


0.35 


2.82 


1.5 


A(eV) 


0.34 


0.80 


0.11 


0.41 


0.43 


0.85 


P(eV • A) 


10.5 


9.69 


8.45 


9.01 


9.6 


8.9 


h{A) 


0.84 


4.7 


0.32 


9.1 


0.3 


1.4 


h{A) 


3.55 


9.6 


1.96 


18 


1.38 


3.87 


5(eV • A^) 


27.6 


108 


4.86 


34.3 


13.6 


41.8 


C 


0.1 


0.5 


0.06 


0.004 


0.4 


0.16 



Table 1: Properties of several direct-gap zinc- blende semicon- 
ductors. Here we denote 6 = 6|i®'^ - the spin-orbit parameter 
defined in [2]. The constant C is defined after (|27|) . 
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